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[^^ ■ D?> and I?7-branes intersecting in 2 + 1 dimensions give rise at low energies to A/" = 4 



supersymmetric Yang-Mills theory coupled to defect fermions in the fundamental repre- 
sentation. This theory undergoes a BKT-type phase transition from a conformal phase 
to one in which the fermions acquire a non-zero mass when the 't Hooft coupling of the 



Hi A/" = 4 SYM exceeds a critical value. To study this transition, we continue the param- 

eters of the model to a regime where a gravitational description is valid. We use it to 
calculate the masses of mesons and the phase diagram as a function of temperature and 
chemical potential. We also comment on the relation of our discussion to the transition 
from the non-abelian Coulomb phase to a confining one believed to occur in QCD at a 
critical number of flavors. 



1. Introduction 

Systems of intersecting branes in string theory have proven useful in studying a va- 
riety of non-trivial dynamical phenomena in quantum field theory which are difficult to 
address using other techniques. Holography provides an efficient tool for analyzing such 
systems in particular regions of the parameter space of brane configurations. For appli- 
cations one is often interested in different values of the parameters, but the gravitational 
description provides a qualitative picture of the dynamics. Quantitative predictions can 
then be obtained by systematically including string (a') and quantum (qs) corrections. 

In this paper, we discuss an example of this paradigm. The system we consider is Af = 
4 SYM with gauge group SU{N), coupled to fermions in the fundamental representation 
of the gauge group, which are localized on a codimension one defect. In this system one 
can study the dynamics of fermions with a continuously tunable interaction, governed by 
the 't Hooft coupling of the Af = 4 theory, A = g'^N. As we will see, below a critical 
value Ac, the theory is conformally invariant in 2 -|- 1 dimensions; above Ac, the conformal 
symmetry is dynamically broken and the fermions acquire a mass. 

This defect QFT can be realized in string theory as the low energy theory on N D3- 
branes intersecting one or more DT-branes along an IR ' . Some aspects of the dynamics 
of this brane configuration were studied in [|I| , and its relation to the Quantum Hall effect 
was discussed in [|]-|] (see also |^-§]). Our main interest in this paper is the physics in the 
vicinity of the phase transition at A = Ac, which does not seem to have been addressed. 

Similar transitions have been discussed extensively in the context of 3 -(- 1 dimensional 
non-abelian gauge theories,a where they are expected to separate the chirally symmetric 
non-abelian Coulomb phase from the chiral symmetry breaking confining one. E.g., in 
QCD with gauge group SU{N) and F fiavors of fermions in the fundamental represen- 
tation, such a transition is expected to occur at a finite value of x = N/F in the limit 
A^, -F — 7- oo. According to this scenario, for x < Xc the infrared theory is conformal and 
preserves the chiral SU{F)l x SU{F)fi global symmetry. For x > Xc the chiral symme- 
try is broken, and the fermions acquire a mass /x, which scales near the transition like 
fx ~ AQCD^wi—d/y/x — Xc). The precise value of Xc is not known but is believed to be 
around 1/4. The large hierarchy between fi and Aqqd is of interest in the context of 
(walking) technicolor (see e.g. [0 for a review). 



See e.g. Q for a recent review and further references. 
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Phase transitions such as the ones described above are sometimes referred to as Con- 
formal Phase Transitions (CPT's) |TI[]. The authors of []T^ pointed out that they occur 
generically when two fixed points of the RG approach each other and annihilate. This 
naturally leads to the above scaling of the order parameter, which we will refer to as Mi- 
ransky scaling or, following [12], as Berezinskii-Kosterlitz-Thouless (BKT) scaling. Some 



recent discussions of holographic models which exhibit such scaling appear in [1T3| - |T7[| 



In the defect QFT that we will study, the dynamically generated mass of the fermions 



is expected to scale like A[/yexp(— 6/-\/A — Ac) as A — > Ac, where h.uv is the UV cutoff and 
h a constant. Since Ac is of order one (as we will see), superficially we cannot study the 
vicinity of the transition using either weak or strong coupling techniques. One approach 
for addressing this problem is to change the parameters such that the transition occurs at 
weak coupling, i.e. Ac ^ 1. Indeed, if one takes the fermions to be localized on a IR ''~'^' in 
IR ' , one finds |T^ that for small e, Ac ~ e^. Of course, one then has to continue the results 



to e = 1 to get quantitative information about the original transition. This can in principle 
be achieved by including higher order perturbative contributions to the dynamics. It gives 
rise to the celebrated e-expansion that plays an important role in the study of second order 
phase transitions in statistical mechanics, and in many other contexts. 

In this paper we will introduce a complementary approach for studying the transition, 
which uses holography. As we will see, one can continue the parameters of the model to a 
region where the transition occurs at a large value of the 't Hooft coupling. There, A^ = 4 
SYM is described by IIB supergravity on AdS^ x S^, and the fermions by probe branes 
propagating in this background. The study of the physics in the vicinity of the transition 
amounts then (at leading order in 1/A^) to analyzing the dynamics of the probe branes 
in AdS^ X S^ . We will use this description to calculate a number of observables, such as 
masses of scalar and vector mesons and the phase diagram of the model as a function of 
temperature and chemical potential. 

In this approach, one does not have to change the dimension of the defect on which 
the fermions are localized. Instead, one changes the number of 3 + 1 dimensional scalar 
fields which couple to the defect fermions via Yukawa interactions. As in the case of 
the e-expansion, one has to include corrections to the gravity picture - in this case, a' 
corrections - in order to study the CPT quantitatively. However, we will see that the 
leading approximation already provides useful qualitative, and perhaps even quantitative, 
information. 



The plan of the paper is as follows. In section 2 we introduce the D3 — D7 brane system 
and its low energy field theory description. We discuss the dynamics at weak and strong 
coupling, and show that this theory undergoes a phase transition at a finite value of the 
coupling, from a conformal phase to one in which the conformal symmetry is dynamically 
broken and a mass is generated for the fermions. 

In section 3 we generalize the discussion to other intersecting brane systems. We 
tune the parameters of these systems such that the phase transition occurs either at very 
weak or very strong coupling, where it can be studied using weak coupling field theory 
techniques and holography, respectively. We show that the transition is conformal in both 



regimes. This can be viewed as due to the mechanism described in |T^ ; as one approaches 
the critical coupling from below, two fixed points of the RG come together, annihilate and 
move off into the complex plane. This "topological" nature protects the transition as one 
varies the parameters, and makes it plausible that it is conformal in the system of section 
2 as well. 

In sections 4, 5 we use the gravitational description to study the system near the tran- 
sition. Section 4 contains an analysis of scalar and vector mesons. As usual in holographic 
systems, we find that in the gravity regime the mesons are very tightly bound, i.e. their 
mass is much lower than that of the constituent fermions. We also find a finite suppression 
of the mass of the lowest lying scalar meson relative to those of other mesons. This scalar 
can be thought of as a pseudo-Goldstone boson of (explicitly) broken conformal symmetry. 
In section 5 we discuss the phase structure of the model as a function of temperature and 
chemical potential. We find that at small T, Hch the symmetry is broken, while at large 
values of these parameters it is restored. The two phases are separated by a line of first 
order phase transitions. 

In section 6 we discuss the conformal phase transition in QCD from our perspec- 
tive. Further comments on our results and possible extensions appear in section 7. Some 
technical details of the analysis are summarized in two appendices. 

2. The D3/D7 system and its low energy dynamics 

The brane system we consider consists of A^ DS-branes intersecting F DT-branes in 
2-1-1 dimensions. We will take N to be large, and F = 1; it is easy to generalize the 
discussion to larger values of F. The branes are taken to be oriented as follows: 

0123456789 
D3 X X X X (2.1) 

D7 XXX xxxxx 



The low energy spectrum of this system includes two sectors. Open strings ending on the 
threebranes give rise to A/" = 4 SYM with gauge group U{N), living in the spacetime labeled 
by (x°,a;^,a;^,a;^). Strings stretched between the threebranes and the sevenbranes give 
charged matter fields localized at x^ = 0. The spectrum of these fields can be determined 
by a worldsheet calculation, but a quick way to find it is to note that the brane system ( |2.1| ) 
is T-dual to the well studied D4/D8 system that plays a role in the Sakai-Sugimoto model 
||18|| . Hence, the only massless fields are 2 + 1 dimensional fermions which transform in the 
fundamental representation of U{N), and can be thought of as chiral (Weyl) fermions in 
3 + 1 dimensions, dimensionally reduced to 2 + 1 dimensions. Strings with both ends on 
the sevenbranes give 7 + 1 dimensional massless fields, which can be treated as classical. 

We will focus on the dynamics of the defect fermions due to their coupling to the 
A/" = 4 SYM fields. In the rest of this section we will discuss the dynamics first at weak 
coupling, where they can be studied using standard field theoretic techniques, and then at 
strong coupling, where one can use holography. 

2.1. Weak coupling 

For small A, the low energy dynamics of the branes is governed by the action 

S = Sn=4 + / d^x{nffi; + g^i^cj)^). (2.2) 

The first part is the action of A/" = 4 SYM, 

j d^xTr Ql)^(/>9p'^(/)9 - \f^uF^'' + ...") , (2.3) 



Sn= 



N=4, 



where cj)^ is a scalar field which transform in the adjoint representation of U{N) and 
parametrizes the fiuctuations of the threebranes in the x'^ direction. 

The second term in ( |2.2| ) is localized on the defect (at x^ = 0), and describes the 
dynamics of the fermions. Note that only the polarizations of the gauge field along the 
intersection couple to the fermions - the field A^ does not couple to them directly. The 
Yukawa interaction is related to the fact that if we displace some or all of the D3-branes 
relative to the DT-brane in the x^ direction, thereby giving some of the eigenvalues of the 
N X N matrix (p^ non-zero vevs, the corresponding components of the fermions get a mass 
proportional to these vevs. This also shows that there are no Yukawa interactions of the 
fermions with the other adjoint scalars (/>*, i = 4, 5, 6, 7, 8. 



The brane configuration ( |2.1| ) preserves an SO{2, 1) x SO{5) subgroup of tlie ten 
dimensional Lorentz group. Tliis symmetry is manifest in the action ( ^.21 ): the first factor 
is the 2 + 1 dimensional Lorentz symmetry of the defect, while the second is the subgroup of 
the SO{6)ii symmetry of A/" = 4 SYM preserved by the defect. The action (|2.2D furthermore 
preserves 2 + 1 dimensional conformal symmetry. In the large A^ limit, this symmetry 
remains unbroken in the quantum theory up to a critical value of the coupling, Ac — 1. 

A common approach to studying the long distance dynamics of fermions in quantum 
field theory is to analyze the infrared properties of their self-energy S(p) (see e.g. ||T9| , pO| ), 
which can be thought of in terms of the non-local gauge invariant observable 



OW{x, y) = ip{x)Vexp 



y 
ig I A-dl 



i^{y). (2.4) 



OW is an Open Wilson Line (OWL) of the sort that plays a role in gauge-gravity duality 
(see e.g. PT[). Its definition depends on the choice of contour going from x to y; a, natural 



choice is the straight line connecting the points x and y [^. The vacuum expectation 
value oiOW can be thought of as a gauge invariant generalization of the fermion two-point 
function that is used to define the self-energy E. As we discuss next, it can be used as an 
order parameter for the conformal phase transition. 

Suppose the dynamics of the defect fermions described by ( |2.2| ) is conformal for A < Ac, 
and massive for A > Ac (we will justify this assumption below). For A < Ac, the one point 
function of OW is then constrained by symmetry to have the forma 

where A(A) is the scaling dimension of OW . In the free theory it is given by A(0) = 2, 
but in general, it is expected to depend on A. 

For A > Ac, the conformal symmetry is broken, and the behavior (|2.5|) is modified. 
The leading correction can be evaluated by applying the operator product expansion to 

lim OW(x,y) ~ i^^— -^^^^ + _____v,v,(.) + . . . (2.6) 

where the ellipsis stands for contributions of higher dimension operators such as ipV^ij), 
and Ab is the scaling dimension of the fermion bilinear operator ^ip^ip. We will write it as 

A,(A) = A(V^V') = 2 + 7(A) (2.7) 



Here and below we omit some overall numerical factors. 



where 7 is the anomalous dimension of i()i(j. 

Taking the vacuum expectation value of ( |2.6D one finds 



In the conformal phase A < Ac, the vev (V'V') vanishes and one recovers ( |2.5| ). For A slightly 
above Ac, this vev is non-zero but very small. Hence one can treat the second term in ( p.8| ) 
as a small perturbation, for sufficiently small \x — y\. 

Another way of breaking the conformal symmetry is to add to the action a mass term 
for tj), m j d^xt/jt/j. For small m, the contribution of this term to ( |2.8| ) can be included 
perturbatively. One finds 

(0W(., y)) ~ '^ (,_^|^(.,-i + (,4t^^W-A. + (x-y)A")-3.A. + ■ ■ ■ ■ (2-9) 

To bring ( |2.9| ) into a more familiar form, it is useful to rewrite it in momentum space. We 
can parametrize the momentum space one point function as 

{OW{p)) ~ p^W-^ (^ + M{p) + ■■■), (2.10) 

where _ 



and /x is a renormalization scale. The one point function (|2.10| ) is closely related to the 
fermion two point function: 

Sf{p) = -, ^T^. (2.12) 

f — m — Ti[p) 

E(p) can be calculated perturbatively by summing IPI Feynman diagrams. It can be 
parametrized by 

E(p) = (1 - A{p))i> + B{p) - m, (2.13) 



which leads to 



with 



^^^' A{p)^-B{p) A{p)p^-M^{py ^ ■ ' 



M(p) = |M. (2.15) 
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At large momenta one can expand ( p.l4 ) as 



(2.16) 



which has a very similar form to ( |2.10|) . A{p) is related to the anomalous dimension of the 
fermion '0 and is typically gauge dependent |jl^,|Ol; its gauge invariant analog in ( |2.10| ) 
is the factor p2~^(-*^). The large momentum behavior of M{p) was discussed in [|l^. An 
analog of the operator product expansion ( |2T^ ) for 'i/'(a;)'i/'(y) leads again to an expansion 
of the form (|2.11|) . Thus, the leading large momentum behavior of M[p) in (|2.16|) is given 
by (|2.11|) and in particular is gauge invariant pO| . 

Note that M(p) ( p.ll|) is a solution of the second order equation 



±{p^M'{p)) + C{X)M{p) 



(2.17) 



where the function C(A) is related to the anomalous dimension 7(A) (|2.7|) , (|2.11|) via 



7(A) = -^ + a/^-C(A), 



(2.18) 



and can be calculated perturbatively in A using standard techniques. The one loop calcu- 
lation described in appendix A yields 

5A 



C{\) 



1271^ 



+ 0(A^ 



(2.19) 



At weak coupling C(A) is small, but as the coupling increases, it can potentially become 
of order one, and in particular approach the value C(Ac) = 1/4 at which the anomalous 
dimension (|2.18|) becomes complex. As we discuss next, this is the location of the CPT. 

In order to solve the second order differential equation (|2.17|) for the self-energy M(p), 
we need to specify the boundary conditions. It turns out that they are 



M'(0) = 0, M(A 



UV) 



0, 



(2.20) 



where h.uy is the UV cutoff of the field theory. The first of these follows from the assump- 
tion that M{p) is smooth near p = 0. The second is the requirement that the bare fermion 
mass term in the Lagrangian vanishes. 

There is a subtlety regarding the boundary conditions (|2.20| ), and more generally 
the definition ( p.lOj) of M. We defined this quantity as the leading deviation from the 
conformal behavior ( p.5|) . Thus, eq. (|2.17| ) is valid only for M{p) ^ p. This condition 
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breaks down at energies on the order of the dynamicaUy generated mass /U, or exphcit 
mass m. In order to impose the first boundary condition in (|2.20| ), we need to extend the 
definition of M{p) beyond the hnear regime described by eq. ( |2.171 ). This is discussed at 
weak couphng in Appendix A, and at strong couphng later in this section. GeneraUy, since 
the solution for M is not expected to change much for < p < /i, one can impose the 
infrared boundary condition at p ~ /x. 

We are now ready to discuss the solution of ( |2.17| ), ( p.20|) . It is easy to see that as 



long as 7 ( 2.18 ) is real (which, assuming C(A) is a monotonically increasing function of A, 
occurs when A < Ac), the only solution of this (linear) equation with the right boundary 
conditions is M = 0. However, for C(A) > 1/4 {i.e. A > Ac), non-trivial solutions exist. 
Defining 

^(A)=C(A)-i, (2.21) 



the general solution of ( |2.17 ) takes the form 



2 



M{p)=An(^\ sin( v^ln-+0j . (2.22) 

The scale /i as well as dimensionless parameters A and 4> ^^^ at this point free. Two of 
them, e.g. A and cp, can be fixed by the boundary conditions (|2.20| ). The renormalization 
scale n can be chosen to be equal to the dynamically generated scale, 

H = M{0). (2.23) 

Any other value is related to this one by a renormalization group transformation. 

Assuming that /i ^ Auv and imposing the second boundary condition in ( p.20|) leads 
to the relation 

v^ln-^ + (/)= (n + l)7r; n = 0,l,2,---. (2.24) 

The solution ( |2.22| ) with n = starts at M = at Auv, and monotonically increases to 
M ^ fj, at p = fi. There the linear approximation M{p) <^ p breaks down, but as mentioned 
above, M remains approximately constant over the remaining momentum interval. The 
dynamically generated scale is given by 

;U~A[/vexpf — ^j. (2.25) 



To have a large hierarchy of scales between fj. and Auy, one must be near the transition, 
where k ^ 1. 



The solutions with n > 1 have a similar structure, but M{p) has in these cases n nodes 
in the interval < p < l^uv- As we will see below, these solutions have higher energies 
then the symmetry breaking vacuum, and are not even locally stable (in the regime we 
will discuss). Hence, they do not seem to play an important role in the dynamics. 

To summarize, as the 't Hooft coupling A increases, C(A) increases and the dimension 
Ab of ?/^'i/' decreases. The CPT occurs when C(Ac) = 1/4. For A < Ac, A^ is real and 
M[p) vanishes; for A > Ac it is complex, and the lowest energy solution has M{p) ^ 0, 
with large momentum behavior given by ( |2.22| ). At A = Ac the anomalous dimension 
( p.l8| ) is 7(Ac) = —1/2, the operator ipip has dimension 3/2, ('0'0)^ becomes marginal,0 
and the two leading terms in the large momentum expansion of M given in ( |2.11|) have 
the same momentum dependence (up to logarithmic corrections). This is the analogue of 
the statement in four dimensional QCD that the transition from the non Abelian Coulomb 
phase to the confining one is expected to occur when /S.[ip%lj) = 2 and the quartic fermion 
coupling becomes marginal [p!9| . 



This picture is also in agreement with the discussion of ||T^ , illustrated in figure 1 . For 
A < Ac, the theory has two fixed points that differ in the value of the four Fermi coupling. 
The anomalous dimensions of il^tl) at the two fixed points takes the values 



7±(A) 



1 

--± 
2 



i - P(^)- 



(2.26) 



As A — 7- Ac the fixed points approach each other, and above Ac they move off the real axis. 
The resulting transition satisfies BKT scaling ( p.25|) . 



coupling 



Fig. 1: Fixed point structure as a function of the 't Hooft coupling A. 



Note that we are assuming here that A.{{iIj%1>)'^) = 2A{7p%p). This assumption is valid in 
the large A^ limit, where our discussion takes place, but in general this relation receives 1/N 
corrections. We will comment on such corrections in the discussion section. 



In the discussion above, we have assumed that the anomalous dimension of the mass 
operator gets to —1/2 at a finite value of the 't Hooft coupling, or in other words that C(A) 
gets to 1/4 at a finite value of A = Ac- The perturbative calculation of appendix A which 
leads to (|2.19 ) is clearly not enough to establish that. However, holography enables us 



to analyze the system, and in particular to calculate C(A), in the strong coupling regime. 
We next turn to that calculation, which will help to determine whether a phase transition 
indeed takes place. 

2,. 2,. Strong coupling 

For large A, the dynamics of A/" = 4 SYM is described by type IIB supergravity (more 
generally IIB string theory) in AdS^ x S^. The L>7-brane ( p.l| ) can be viewed as a probe 



propagating in this background [^; its back- reaction on the geometry can be neglected 
to leading order in 1/A^. To describe the intersecting brane system ( p.l|) it is convenient 
to write the metric on AdS^ x S^ as follows: 

ds^ = (jY dx^,dx^ + (-^ {dp^ + p'^dnl + {dx^f), (2.27) 

where n = 0,1,2,3, and L is the AdS radius,a L^ = 27r^A. The six dimensional space 
transverse to the threebranes has been split into an IR corresponding to the directions 
(45678), which is described by the spherical coordinates (p, O4), and x^. Thus, these 
coordinates satisfy p^ = (x^)^ H — ■ + (x^)^, r^ = p^ + (x^)^. The background is supported 
by RR five-form fiux, but this will not play a role below and we will not write it explicitly. 
In the brane configuration (|2.1|) , the L>7- brane is stretched along the surface x'^ = 
x^ = 0. The induced metric on its worldvolume is given by 

ds^ = f^YdXadx" + (-^ {dp^ + p'^dQl) , (2.28) 



where a = 0, 1, 2 are directions along the D3/D7 intersection. The metric (|2.28|) is that of 
AdS4 X iS^ ; it makes manifest the conformal symmetry of the problem, as well as the other 
symmetries mentioned above. However, this configuration is unstable to condensation of 
the scalar field corresponding to the fluctuations of the L>7- brane in x^. Indeed, if we 



^ We set a' = 1. 
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generalize the ansatz for the shape of the brane from x^ = to x^ = /(p), we find the 
induced metric 

ds^ = (^^y dxadx'' + (^\ {[l + f'{pf]dp'' + p^dnl). (2.29) 

The profile f{p) is determined by minimizing the DBI action 



So7 = jd'xj dp-^^j^pWl + fip? (2.30) 

where the prime denotes differentiation w.r.t. p and we omitted an overall multiplicative 
constant. As is usual in holography, one should think of the radial direction p as corre- 
sponding to energy scale in the field theory. The upper limit of the integral in ( |2.3C1| ), A, 



is related to the UV cutoff in the field theory via [p3| , p^ : 

A ~ Auv^- (2.31) 

The Euler-Lagrange equation that follows from ( p.29|) is 

d ( P^ f\p) \ . '2f{p) 4 /T— 77732 n to 'io^ 

The ground state of the system corresponds to the lowest energy solution of this equation 
which satisfies the boundary conditions 

/'(O) = 0, /(A) = 0. (2.33) 

The first is necessary for regularity of the shape of the brane near the origin of the IR 
labeled by (45678); the second is the statement that we do not displace the sevenbranes 
from the threebranes by hand, i.e. that the bare fermion mass vanishes. 

The simplest solution of (|2.32|) with these boundary conditions is f{p) = 0. It preserves 



the scaling symmetry p — )■ ap, / — )■ af of (|2.32|) , and corresponds to the conformally 
invariant state (|2.28 ). However, as we will see next, the ground state of this system has 



f{p) 7^ 0, and breaks the conformal symmetry. 

To find it, we start at large p, where we expect f{p) to be small. Expanding ( p.32| ) in 
f{p) yields the linear equation of motion 

^(pV'(p))+2/(p)=0. (2.34) 
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If we think of f{p,x"') as a scalar field in AdS^, ( |2.34| ) is essentially the Klein-Gordon 
equation it satisfies. It is easy to check that the mass of this field is below the Breitenlohner- 
Freedman (BF) bound [^. Thus, the system is unstable to its condensation^ 



The general solution of (|2.34|) is 




f{p) = A^i|^^] sin(^ln^ + <^|, (2.35) 

where we introduced two dimensionless parameters A and (p, and an arbitrary scale p. 
Solutions with different values of /i are related by the scaling symmetry mentioned above. 
The result (|2.35D is reliable at large p; as p decreases, the small f{p) expansion eventually 



fails, and we have to go back to the full equation ( p.32|) . As in the previous subsection, we 
choose the scale p to be 

p = /(O) (2.36) 

and then fix the parameters A and (p by imposing the boundary condition /'(O) = 0. From 
( p.36| ) we find (numerically) A ^ 0.761; (^ ~ 1.086. To fix /(O), we need to impose the 
second boundary condition in ( |2.33| ), /(A) = 0. This gives an infinite number of solutions, 

//^ = AexpU^((7i + l)7r-(/))j, n = 0,l,2,--- (2.37) 

some of which are shown in figure 2. The nth solution has n nodes in the interval < p < A. 
To identify the ground state of the system, we need to compare the energy densities 
of the solutions that we found, which are the values of the action ( |2.3U| ) evaluated on the 
solutions. The energies of all the states behave at large A like L^A"^/3, but the differences 
of energies are finite as A — )■ oo. One finds that the energy grows with n, so the solution 
with n = (no nodes) is the ground state. The conformally invariant solution /(p) = 
can be thought of as the limit as n — t- oo of the others; it has the highest energy. If we set 
the energy of this solution to zero, the ground state has energy density 



'^0 



-Q.np^L^/gs ~ -p^N/VX. (2.38) 
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See e.g. |26-^| for other discussions of tachyons below the BF bound in AdS. In these papers 



the infrared instabilities due to the presence of these tachyons were associated in the dual field 
theory to double trace operators. In our case the tachyon is a single trace operator corresponding 
to a light mode of an open string whose ends lie on a iD7-brane in AdSs x S and the role of the 
double trace operators is played by (V'V')^- 
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Fig. 2: Solutions of ( |2.32| ) with A = 4.7305 and fi = 1. The solid, dashed and 
dotted lines correspond to n = 0, 1,2. The n = 1,2 solutions have been rescaled 
by factors of 8 and 64, respectively. 



To compute the energy density in 



we restored the factor of I/qs in the normahzation 



of the DBI action ( |2.30| ) . The A^ dependence of the vacuum energy is compatible with the 
fact that we are dealing with a system of A^ fermions, and the symmetry breaking affects 
their state. The A dependence suggests that another scale, in addition to fx, plays a role in 
the dynamics. As we will see later, the two important physical scales are ^ and n/yX. The 
former corresponds to the dynamically generated mass of the fermion, given by the energy 
of an open string which stretches at p = between the DT-brane (located at x'^ = n) and 
the horizon at x^ =0. The latter is associated with the meson masses and the temperature 
of the phase transition to the symmetric phase. 

From the brane perspective, one can think of the energetics as due to a repulsive force 
between the D3 and 1^7-branes, which (at strong coupling) pushes the sevenbranes away 
from the origin in x^. This breaks the x^ — ;■ —x^ parity symmetry, as well as conformal 
symmetry, and generates a mass for the fermions. The analysis of the previous subsection 
shows that for small A this repulsive force is too weak to generate a mass. 

Our treatments of the fermion self-energy M{p) in the previous subsection, and of the 
sevenbrane profile f{p) in this one, were very similar. Indeed, the linear equation ( |2.17| ) 
is a direct analog of ( |2.34| ), with M ^ f and p ^ p. The boundary conditions ( |2.2C1| ) and 
( p.33| ) are the same as well. The large momentum expansion of M ( p.ll| ) is mapped to 
the standard large p expansion of a scalar field in anti de-Sitter space, with the leading, 
non-normalizable, term corresponding to a coupling in the Lagrangian of the dual field 
theory, and the first subleading correction encoding the vev of the corresponding field. 
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The boundary field t/jtl^ whose couphng and expectation value appear in ( 2.11 ) is 
related by holography to the bulk field /(p). Hence it is natural to think of the large p 
expansion of f{p) as providing information about the (strong coupling limit of the) large p 
expansion of M{p). In particular, comparing ( p.22|) to (|2.35|) we see that limA^oo ^(A) = 
7/4, i.e. C{X —7- oo) = 2. Assuming that «; is a continuous function of A, it has to pass 
through zero at a critical coupling Ac- 

To summarize, we see that the system ( |2.1| ), (|2.2| ) undergoes a phase transition at a 
finite coupling, which is out of reach of both the weak and strong coupling expansions. 
Indeed, for small A we found M{p) = and no dynamically generated scale, whereas for 
large A the dynamically generated scale, po in (|2.37| ), is comparable to the UV cutoff A. 
The transition could either be continuous, with the dynamically generated scale vanishing 
for A < Ac (defined by k(Ac) = 0), and continuously approaching zero as A —t- Ac from 
above, or it could be discontinuous, with the scale in the broken phase bounded from below 
by a value of order the UV cutoff. In the former case the transition must be conformal, 
driven by the mechanism described in the previous subsection, and exhibit Miransky scaling 
( pj.25| ). In the latter, the transition must occur at a coupling strictly below Ac- 

In order to distinguish between the two possibilities, we will describe in the next 
section two perturbative expansions that allow one to analyze the physics near the phase 
transition. By changing the parameters of the system ( |2.2| ), we will be able to push the 
transition towards weak or strong coupling, and show that in both limits the dynamically 
generated scale can be made arbitrarily small in its vicinity. This provides strong evidence 
that the phase transition in this system is conformal. 

3. Perturbative expansions 

In the previous section we saw that the D3/D7 system undergoes a phase transition 
at a coupling of order one. In this section we will discuss approximation schemes that 
allow one to study this transition. To this end, we will replace the flavor D7-brane by a 
Dp-hrane which shares d spacetime dimensions with the -D3-branes, and is extended in 
n = p + 1 — d additional spatial directions. We take the orientation of the branes to be 

0123456789 
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where solid (dashed) hnes indicate directions along (transverse to) the branes. The original 
D?,/D7 system ( pTT| ) corresponds to {d,n) = (3,5). The supersymmetric D3/D5 system 
studied in [0,|31|] corresponds to d = n = 3, and one can consider many other systems 
with various d and n. We will eventually take these parameters to have general non-integer 
values. For d, this is familiar from the e-expansion, but we will find it useful to vary n as 
well. 

The brane configuration ( p.l|) is invariant under 0{d—l, 1) xO{4—d) xO{n) xO(6 — n). 
In addition to the A/" = 4 SYM fields living on the -D3- branes, the low energy spectrum 
contains the massless states of open strings stretched between D3 and Dp-hranes. The 
number of directions of space in which these strings have Dirichlet boundary conditions on 
one end and Neumann on the other is i^ND = 4 + n — d. If this number is smaller than 
four {i.e. for n < d), the lowest lying (bosonic) excitation in this sector is tachyonic. The 
dynamics of open string tachyons is well understood; their presence typically means that 
the brane configuration can continuously decrease its energy by classical rearrangement 
of the branes. Our interest here is in different dynamical phenomena; therefore we will 
consider only n > d where such tachyons are absent. 

For n = d, the brane system ( |3.1|) is supersymmetric; it preserves eight of the 32 
supercharges of type IIB string theory. The massless spectrum contains in this case a 
hypermultiplet which transforms in the fundamental representation of SU{N). An example 



of this is the system considered in [30,31|. Since our interest here is in the dynamics of 



defect fermions, we will restrict to the case 

n> d (3.2) 

in which the lowest lying bosons in the (3,p) sector are massive, with m^ ^ n — ci (in string 
units). The fermions can be thought of as chiral spinors of SO{d + 5 — ti, 1) decomposed 
under SO{d — l, 1) x 5*0(6 — n); they are singlets of 5*0(4 — d)x SO{n). We will eventually 
be interested in continuing d, n to general non-integer values, but to perform calculations 
at weak coupling it is convenient to take d — n E 2Z and continue the results. 

Our goal in the rest of this section is to repeat the discussion of the previous section 
for all d and n and look for regimes in which the critical coupling Ac is very small or very 
large, so we can use perturbation theory or holography. 
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3.1. Weak coupling 



The low energy effective field theory corresponding to the brane configuration (|3.1 
described by the action 



S = Sn=4 + f d'^xiiiPfiJj + gi;Tm(p"'iP) , (3.3) 



where the vector index m runs over the 6 — n transverse directions in ( p.l| ) . (f)"^ are four 
dimensional fields that transform in the adjoint of SU{N) and parametrize fluctuations of 
the threebranes in these directions. They are generalizations of the field (p^ in the previous 



section, and the origin of the Yukawa couplings in (|3.3| ) is the same as there. 

The weak coupling calculation of the fermion self-energy ( |2.12| ) - (|2.16| ) is described 
in appendix A. One finds that eq. ( |2.17| ) generalizes to 

^ {p^-'M'{p)) + C{X)p^-'M{p) = 0, (3.4) 

where 

C{X) = j^{2{d - 1) + (6 - n)) + 0{\^). (3.5) 

The solution of ( |3^ ) at finite coupling depends on the sign of 

k{\) = C{\)-{^^ . (3.6) 

For K < {i.e. for sufficiently weak coupling), it takes the form 

where 



d-2 Ud-2^^ 



7(A)=7+(A) = _ + W^_ -C(A) (3. 



is the anomalous dimension of t/jt/j. In the limit m — )■ the theory becomes conformal, and 
M(p) vanishes, as in the discussion of section 2. 

In this range of couplings, we also expect a second fixed point, which is formally 
obtained by adding to the Lagrangian the irrelevant double trace operator {"ipipY ^^d 
fiowing up the resulting RG trajectory (see figure 1). At this fixed point the anomalous 
dimension 7 takes the form 



d-2 lfd-2^^ 



7(A)=7-(A) = —-\j[^^] -C{\). (3.9) 
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Hence, the operator (^^)^ is relevant (for d > 2); adding it to the Lagrangian leads 



back to the fixed point (|3.3[ ). Note that the dimension of i/ji/j at the new fixed point, 
Auyi^i/ji/;) = (i— 1 +7_, satisfies the unitarity bound A > {d — 2)/2 for all d < 4 (assuming 
C(A) > 0). This is consistent with the existence of such a fixed point for all values of the 
coupling. 

For A = this fixed point is familiar from studies of the NJL model. As is well known 



(see e.g. [^2| for a review), in d < 4 dimensions this model is solvable in the large A^ limit. 
Using standard large N techniques one can show that the fermion bilinear a = i/ji/j has 
dimension one in the UV, in agreement with ( |3.9| ) (with C(A = 0) =0). For A 7^ the 
model is no longer solvable, but it is natural to expect that the NJL UV fixed point can 
be smoothly continued to this regime. 

When the coupling exceeds a critical value Ac for which C(Ac) = {d — 2)^/4, the 
two fixed points discussed above merge and move off the real axis []T^. The anomalous 
dimensions (|3.8D , (|3.9[) become complex, k turns positive, and the solution of ( |3.4|) assumes 



the form 

M{p) = An(-j ' sm(^/^ln- + <p] . (3.10) 

In this regime, conformal symmetry is dynamically broken. A and (p can be calculated as 
in section 2, and the scale fi can be set equal to the dynamically generated scale, as in 
(PI). 

The transition between the two regimes occurs at the point where k (|3.6D vanishes. 
In section 2 we discussed the case (i = 3, where this happens at a coupling of order 
one. However, we see from ( |3.5| ), ( |3.6[ ) that if the dimension d is close to two {d = 
2 + e), the critical coupling can be made small (Ac ~ e^). In this case, the weak coupling 
approximation is valid in the vicinity of the transition. 

It is easy to understand why near two dimensions the CPT happens at weak coupling. 
In general, the transition occurs when the anomalous dimension ( p.8|) takes the value 
7(Ac) = —{d — 2)/2, so A{'i/j'i/j) = (i/2, and ('(/''i/')^ is marginal. As is well known, for d = 2 
this operator is marginal already in free field theory. Thus, for d slightly above two, the 
anomalous dimension needed to make it marginal is small. 

To summarize, we see that the CPT of section 2 can be studied using weak coupling 
techniques by continuing the dimension d to 2 + e and using the e expansion [jl2|. Of 
course, to obtain quantitative results in d = 3 one needs to continue to e = 1, but already 
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in leading order in e one finds a phase transition which is quahtatively similar to what one 
expects; thus, the higher orders in e are not expected to change the qualitative picture. 

For < e, K <^ 1, there is a large range of momenta /x ^ p ^ A in which M{p) is well 
approximated by the solution of ( p.4]) with k = {i.e. C{X) = {d — 2)^/4). In this case 
the two terms in ( |3.7[ ) have the same momentum dependence, and the general solution of 
has the form 



d-2 



^(P) = ^(^) ' fCiln^ + CsV (3.11) 

where Ci, C2 are dimensionless constants. From the point of view of ( |3.10| ), this is a regime 
in which the argument of the sine is small. Comparing coefficients, we conclude that for 
small K 

A = ^, <P=^V^. (3.12) 

^/K Ci 

The dynamically generated scale can be obtained by imposing the boundary condition 
M(A) = ( 12:201 ). Looking back at ( pl| ), one finds the Miransky (or BKT) scaling 
behavior 

M(0) = ;U ^ Aexp ( — ^ ) . (3.13) 



Note that the fermion mass m^ is not exactly equal to /U, but is comparable to it. Indeed, 
m^ is obtained by solving the equation M(to^) = m^ (see (|2.14|) ). The self energy M{p) 



has the following structure: it is roughly constant and equal to fi for p < fi, and then 



crosses over to the behavior (|3.11| ) for p > fj,. p = fx is in the crossover region between the 
two behaviors, but since both give M{fj,) ~ fj, there, it must be that m^ ^ /i. One can 
obtain the ratio m^/fj, more precisely by solving the equations of appendix A numerically. 

3.2. Strong coupling 

To extend the discussion of subsection 2.2 to general d, n, it is convenient to write 
the metric of AdS^ x S^ as 

ds"" = (^^y dx^dx^" + (^^ {dp^ + p^dnl_^+df + fdnl_J, (3.14) 

where fx = 0,1,2,3, while (p, O^i-i) and (/, Os-^) are spherical coordinates on the IR" 
and IR ~" transverse to the L>3-branes in (|3.1|) . The radial coordinate of AdSs is given by 
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As before, we take the Dp-hrane to wrap IR ' x S"^ ^ and form the curve / = /(p). 
The induced metric on the Dp brane is then given by 

dslp=(^^ydxadx''+(^] [{l + f'{pf)dp^ + p^dnl_,] , (3.15) 

where a = 0, 1, 2, ■ ■ ■ , d—1 runs over the intersection. For / = 0, ( |3.15p describes AdSd+i x 
S'^~^. As in the discussion of the D3/D7 system in section 2, this configuration is in general 
unstable to condensation of /. To study this instability, we write the DBI action for /, 



d — n 



Sdp = jd'^xldpf^^^j^) ' p-iyiTT^. (3.16) 

The equation of motion for / (the analog of (|2.32|) for general d and n) is 



§-(ip' + /')^'"'^^/'-7J= I + (^ - d){p' + fr—'^''fp-'V^Vr = 0. (3.17) 



The boundary conditions are given by (|2.33| ). The conformal invariance of the system 



leads to a scaling symmetry of ( p.l7|) , p -^ ap, f -^ af- As in the D3/D7 system, this 
symmetry is in general dynamically broken. 

For large p, we expect / to be small and slowly varying. In this regime, ( |3.17] ) reduces 
to 

^(p'-V')+(^-^)p'"'/ = 0. (3.18) 



This equation is identical to ( |3.4] ) with the map M^/,p^p, C ^ n — d. The relation 
between M and / was discussed above (in subsection 2.2). Following that discussion, we 
see that the strong coupling analysis implies that 



and {^) 



Coo = lim C(A) = n-d (3.19) 

A— >oo 



d-2^ ^ 



Koo = lim k(A) = n — d — — - — . (3.20) 

A— >-cxD y 2 y 

When Koo < 0, the system remainaS in the conformal phase for all A. The anomalous 
dimension of '^'0, 7(A) (|37^), can be calculated at small A from ( p.5| ), and for large A from 
( P.19|) ; it remains real for all A. 



Assuming that k{X) is a monotonic function. 

19 



On the other hand, when Kqo > 0, i.e. foiQ 

n>nc = d+ (^^^ 1 , (3.21) 



d-2^ ^ 



the system undergoes a CPT at a couphng Ac for which k vanishes. In d — 3 dimensions, 

nc = 3 + -. (3.22) 



As a check, the D3/D7 system (for which n = 5) is indeed in the regime (|3.21D, in 



agreement with the discussion of section 2. Supersymmetric systems, which have n — d, 
are always in the conformal phase Kqo < ( |3.2C1| ) as one would expect. 

In the regime ( p.21| ), the transition from conformal to massive behavior generally 



occurs at a coupling of order one. In order to study the transition, we vary n to 

n = nc + S (3.23) 



such that the asymptotic value of k, Kqo =5^1 ( p.20|) , takes a small positive value 



and analyze the resulting theory in the DBI approximation. We have thus pushed the 
transition into the regime where we can use holography. Of course, this analysis takes 
place in a theory with the wrong value of n; the theory we are interested in has n = 5, 
while the DBI approximation is reliable for n near the critical value ( |3.22|) . In order to 
obtain quantitative predictions, we need to include higher order corrections in 6. These 
corrections correspond to a' corrections to the DBI action (|3.16|) . 



The spirit of the approximation here is similar to that of the e-expansion. There, 
e.g. in order to compute the critical exponents of the three dimensional Ising model which 
are governed by the IR fixed point of (p'^ field theory in three spacetime dimensions, one 
continues the dimension of spacetime to 4 — e, so that the fixed point becomes perturbative. 
In our case, we push the interesting dynamics to strong coupling, where it can be analyzed 
using holography. This leads to a kind of gravitational e-expansion. As in the original 
e-expansion, the hope is that the physics for a particular value of k depends smoothly on 
71, so that the leading (in this case DBI) approximation provides a good qualitative guide 
to the dynamics. 



Note that Uc is always in the range (3.2). 
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To analyze the theory in the regime where Kqo is smaU and positive {i.e. n is only 
slightly above Uc), we need to solve the DBI equation of motion (|3.17|) for /. At large p, 
this equation reduces to ( |3.18| ), whose solution behaves like 

f{p) = A^f^] ' smf./^ln^ + A. (3.24) 

The dynamically generated scale is given by 



TV 



rir- 



/(0) = ^^Aexp ^ , (3.25) 



as in (|3.13|) . As before, the value of the radial cutoff A in the AdS space is related to the 
UV cutoff in the field theory via ( p.31| ). To focus on the physics near the transition, one 
can take the double scaling limit A — t- oo, Koo — j- (i.e. n — ;■ Uc), with fx held fixed. We 
will refer to this as the BKT limit. In this limit, ( p.24| ) behaves for p ^ fi as (compare to 

d-2 

The constants Ci and C2 can be obtained numerically. For the case of interest here, d = 3, 
they are given by 

Ci ~ 0.53575; C2 ^ 0.90421. (3.27) 

As in the weak coupling analysis, one can use them to calculate A and (p in ( |3.24| ), via 

(EH)- 

Note that eq. (|3.26| ) can be thought of as the asymptotic form of a scalar field at the 
BF bound. The leading term (proportional to Ci) can be thought of in the boundary theory 
as due to adding the dual operator, i/ji/j, to the Lagrangian; it encodes the dynamically 
generated mass of "i/^. C2 is dual to the vev {i/jiIj). We see that the breaking of conformal 
symmetry at the CPT is explicit, in agreement with the discussions in the technicolor 
literature and |T^. 



So far we discussed the dynamical generation of a mass in the model (|3.3|) . One 
could also add an explicit mass term to the Lagrangian, which corresponds in the brane 
language to separating the threebranes and sevenbranes in the transverse IR ~'^. At strong 
coupling, this is attained by changing the UV boundary conditions from /(A) = to 
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/(A) = fuv > 0. Let us define fm = /(O) when /(A) = fuv, and /o = /(O) when 
/(A) = 0. The UV boundary condition leads to 



d-2 



fuv = AU[^j ' gin(^y^lnAA + <^). (3.28) 



^ m 



A 



Using the fact that y^K^ln ^ + (^ = tt, we find in the BKT limit 



d-2 

fuv Ci(fy-) ' In 4^, (3.29) 



fm '\A J /o 

where Ci is the constant in ( |3.12| ), ( |3.27| ). The appearance of the UV cutoff A in ( |3.29| ) 
can be understood as follows. We can rewrite this equation as 

fuvA"^ =Ci if r^)hn^. (3.30) 

Jo 

The cutoff dependence in ( |3.3U| ) implies that the coupling fuv^ which naively has dimension 
one, has in fact dimension d/2 in the quantum theory. The operator that couples to fuv 
is the mass operator t/jt/j, and its dimension near the BKT transition is d/2, as discussed 
above. Thus, fuvi^i^ has dimension d, as necessary. 

From the point of view of the holographic bulk description, in the BKT regime the 
scalar /(p) which describes radial fluctuations of the Dp-hrane in the transverse directions 
has a mass slightly below the BF bound. On the other side of the transition, where n 
is slightly below Uc ( |3.21| ), the mass of / is slightly above the BF bound and the AdS 
vacuum (|3.15 ) with / = is stable. As is well known [^|, scalars whose mass is slightly 



above the BF bound can be quantized in two different ways, which lead to two distinct 
CFT's differing by double trace operators. These CFT's were discussed in the previous 
subsection; they give rise to the anomalous dimensions (|3.(j|), (|3.9|). 

Since for Koo > the dynamics generates a mass for the fermions, one expects the 
trace of the stress tensor to be non-zero. To calculate it, we need to evaluate the action 
( p.l6|) on the solution of the equation of motion (|3.17| ). As mentioned in section 2, this 



action is UV divergent, but differences of actions for different solutions are finite. Setting 
to zero the energy of the conformally invariant solution /(p) = 0, one findsQ 



—L^'-'^u'^r^-NXi'^y-^i 



Sq L"-V ~ -NX^ — ) -'n", (3.31) 

9s 



The numerical computation of the energy density in the BKT hmit at strong coupling is 
outlined in Appendix B. 
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where we used the fact that we are in the BKT regime n ^ ric ( p.21|) , and restored the 
factor of I/qs in the normahzation of the DBI action ( |3.16D . Eq. ( |3.31| ) imphes that the 
stress-tensor has a non-zero vev: 

(Tab) = SoTjah- (3.32) 

This is a very reasonable result; in the BKT limit, the vacuum energy can be expressed 
in terms of quantities that remain finite in this limit, the scale ^ and 't Hooft coupling 
A. This result also suggests that the theory contains extra scales, in addition to /x. In the 
next section we will analyze the spectrum of mesons and see that this is indeed the case. 
Finally, we would like to comment on the solutions with nodes, which were briefly 
discussed in section 2 (around eq. ( |2.37[ )). For small k, the solution with n nodes satisfies 



/(O) 



Uri — Aexp I -—^ j ~ (Uexp I 1 . (3.33) 



In the BKT limit, all ^^ with n > 1 go to zero, so all these solutions approach the flat 
solution f{p) = 0. It is also easy to calculate their energy densities. This can be done by 
splitting the p integral in ( |3.1ti| ) into two parts: from p = to the flrst node, and from 
the first node to p = A. The first integral is identical to the one for the ground state, and 
gives 

Sn ~ -NXi^T-^i^t (3.34) 

The second can be shown to vanish using a combination of two facts. One is that beyond 



the first node, the linear approximation (|3.18|) is excellent for small k, and one can replace 



the full action ( |3.16|) by a quadratic one. The second is that on-shell the integral of the 
quadratic Lagrangian reduces to boundary terms, which vanish since f{p) = at the 
boundaries (both of which are nodes). Thus, we see that the actions of the solutions with 
nodes ( |3.34|) all vanish in the BKT limit, in agreement with the fact that they all approach 
the trivial solution. 

4. Spectrum of fluctuations 

In the previous section, we determined the shape of the probe brane f{p) in the BKT 
regime and found that the theory dynamically generates a scale p = /(O). In this section, 
we will study small excitations about this vacuum, which correspond to mesons. We will 
see that the meson masses are set by a new scale p, defined as 
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This relation between the meson mass scale and the fermion mass scale is similar to that 



of other holographic models (see e.g. [34]). 

We first consider cr-mesons, which correspond to fluctuations of the order parameter 
for the breaking of conformal symmetry. The lowest of these is the would-be Goldstone 
boson of broken scale invariance, the analog of the techni-dilaton in walking technicolor 
theories. We will see that it is lighter than the other states in the spectrum, but not 
parametrically so. We will also discuss vector mesons, the analogs of the (techni-) p-meson 
in QCD (technicolor). 

Since the vacuum with /(p) 7^ spontaneously breaks the global 0(6 — n) symmetry 
to 0(5 — n), when n < 5 one expects to find massless Goldstone bosons that parametrize 
the coset 0(6 — n)/0(5 — n). We will describe them and calculate their mass when the 
symmetry is broken explicitly. 

4.1. a-mesons 

To study radial excitations of the probe brane, we expand the radial scalar / ( |3.14| ) 
around the solution of ( p.l7|) , 

/(p,x'^) = /(p) + y(p,x"). (4.2) 

Since the DBI action preserves d dimensional Lorentz symmetry, it suffices to take the 
perturbation y to be a function of p and t. For such configurations, the action takes the 
form 

Sdj> = -J d^x j dp-^^ Vr4-^V + (/' + 2/0^^4. (4.3) 

If 2/ = 0, this agrees with (|3.16|) (up to an overall constant, which we do not keep track of 



here and below). 

Expanding (^]^) around the solution /(p) and keeping only terms quadratic in j/, we 
find the action 



S2 = - j d'x j dp {A{p)y^ + B{p)y'^ - L^C{p)f) , (4.4) 



where 



d-n p--\{d-n-l)f + p^) r-rjr^^ d {n - d)p^-'ff{p + p^' 



B{p) 



2 (/2 + ^2)2--i. dp\ v^TTT^ 



d — n 
I ~2~ 



2(l + /'2)3/2 



C{p) ='- ^-1 ^1' > , 

2v/rT7^ 



(4.5) 
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To find the spectrum, we take y{p,t) to have the form y{p,t) = 'i/'(p)e*"^*. The linear 
equation of motion for '0 that follows from ( [4 .41) is 



— {B{p)^'{p)) + m'L^C{p)i;{p) - A{p)^l:{p) = (4.6) 

with the boundary conditions ( |2.33| ), i^'iQ) = V'(^) = 0- Since (|4.6| ) is a Unear equation, the 
overall scale of the solution is arbitrary. Thus, ([4.6|) is over-constrained and one expects 
solutions to exist for a discrete set of masses m. 

We now specialize to the physical case d = 3. For large p, the equation for the 
fluctuations ( [4.6|) becomes identical to the equation of motion ( |3.18|) for the embedding 



function f{p) in the linear regime. Hence, its solution takes the form (|3.24| ), 

ifjip) ^ A^ p l-j sin f y^ log - + (/)<, j . (4.7) 



In order for '(/'(p) to vanish at the UV cutoff, the phase (pa in ( [4.7D must be exactly equal 
to (f) in (|X^ . 



In the BKT limit, (|4.7|) reduces to (see the discussion around (|3.26D ) 



^(P)^^(^) fcf'^^ln^+C^-^'l. (4.8) 



The values of C^ and C2 can be obtained by numerically integrating eq. ( [4.6|) starting 
from p = 0, with the initial conditions '0(0) = 1; '(/''(O) = 0. Varying the value of m^ leads 
to different values of C}'^ and C2 ; the spectrum of fluctuations is obtained by demanding 
(^^ = (^or (see(|T|), i^m\)) 

Ci''^/C[^^ = C2/C1 ~ 1.6877. (4.9) 



In fig. 3 we plot the results of the numerical integration. Note that eq. ([4.9| ) is never 
satisfied for negative m'^, which is consistent with the stability of the vacuum. The lightest 
few states have 

m'^/p'^ ^ 0.44, 9.65, 26.63, 51.35, 84, ■ ■ • , (4.10) 



where p is definied in ( |4.1| ). This spectrum is well described by the general formula 

m^/p2 ^ 3.89?i2 + 5.3271 + 0.44; n = 0, 1, 2, ■ ■ ■ . (4.11) 
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Fig. 3: &^^ jC^"^ as a function oirn^ jji^. 

It exhibits the characteristic behavior m„ ~ n at large excitation number n. This is similar 
to other holographic models, although here the behavior does not have an obvious Kaluza- 
Klein interpretation. We also see that the typical mass scale of the mesons is 7^ ~ ^/vA. 



Thus, at large A the mesons are deeply bound, as in [p4 |. 

Another notable fact is that the lowest lying meson, which corresponds to n = in 
( [4.11| ), is quite light relative to the others. Its mass is mo — 0.66/1, while the next lightest 
meson has mass ra ^ 3.I7I, and the asymptotic separation between subsequent masses 
in ( [4.11|) is m„,+i — rrin — 1.9771. Thus, mo is smaller by a factor of three to five than 
the typical mass scales in the problem, but the separation between the two scales is not 
parametric in the BKT limit. The lowest lying cr-meson can be thought of as a pseudo-NG 
boson of broken conformal symmetry. Since the breaking is explicit, it is massive, but is 
nevertheless relatively light. 

A similar issue has been discussed in the technicolor literature, where the question is 
whether as one approaches a CPT by tuning the number of flavors in QCD, the a meson 



is parametrically lighter than all the other mesons; see e.g. p5| , |36[] for recent discussions 
and p7| - |39[] for earlier work. There, the answer to this question is unclear, as it requires 
finding the spectrum of light states in a strongly coupled field theory. 

Since many of the arguments for a light dilaton, such as those of ||3^], apply to our 
case, it is interesting to see what happens here. Our results do not support the point of 
view advocated in [^5| that as one approaches the CPT, the ratio of mass of the dilaton 
to that of other mesons goes to zero. However, we did find that the mass of the lightest 
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cr-meson is somewhat suppressed. Interestingly, a similar possibility was suggested in QCD 
in 



pql and references therein. In a certain (uncontrolled) approximation it was found that 
the techni-dilaton is lighter by about a factor of three than the lightest vector mesons. 
We will discuss vector mesons in our model in the next subsection, and will find a similar 
ratio. 

So far we have discussed the radial excitations of the ground state, which from the 
point of view of the system with small but finite k (and thus finite UV cutoff A) is the 
configuration with no nodes in the interval p G (0,A). One can apply the same tools to 
analyze the solutions with nodes discussed in the previous sections. For the solution with 
n nodes, one finds n tachyons with masses of order 

2 -2 „„ f 27r(n 



' f^r 



?7ij ~ — /U exp I — 1 ; i = 1, ■ ■ ■ ,n. (4-12) 



These tachyons signal perturbative instabilities, due to the fact that the system can con- 
tinuously lower its energy by eliminating the nodes. The time scales for this process 
correspond to the masses ( [4.12|) . The fastest time scale, corresponding to z = n in ( [4.12|) , 



is of order /x, as one would expect. The time scales corresponding to i < n are much longer, 
but involve fine-tuned dynamics. 

4-2. Vector mesons 

To compute the spectrum of vector mesons, we turn on the gauge field on the world- 
volume of the probe brane. The DBI action then reads 

S = - f d'^x f dp^/-det{g + 27:F), (4.13) 

where g is the induced metric on the brane. Expanding in the field strength F and keeping 
only quadratic terms, we find 



d — n 



d^ j„ y ^ J \ „n-l 



2 I /2 \ — 



p' + f 



S = So- d'^x dpi ^—rr^ p"" Vl + /' 



F F°'P T^F uF"-^ 

-'- ap-'- i„ 



L2(l + /'2) 4(p2 + / 



2^2 



(4.14) 



where the indices along the intersection a, h are raised with the fiat metric rj"'^ . 

For F = 1 (one fiavor brane) one can choose the gauge Ap = 0. For larger F, the 
dynamics in general breaks the U{F) global symmetry to Yli U{Fi) with ^^ Fi = F, and 
the off-diagonal components of Ap play a role in the analysis of NG bosons. 
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The remaining components of the gauge field can be expanded in an orthonormal basis 

A„(x^p) = 5^S(")(x^)^.(p). (4.15) 

n 

Suppose that the "(pn are solutions of 

-G-^dp{KdpiPn) = mliP^, (4.16) 

where 

n-l 

Kip) = ^ „. 

(4.17) 

„, n-d + 4 ' 
M j 2 

and we impose the orthonormality condition 

/ dpG{p)ip^{p)ipn{p) = Smn- (4.18) 

Then integration by parts leads to 

/ dpKdpl/jrndpiJn = mlSrnn (4-19) 



K{p) = 
G{p) = 


pu-l 


{p' + p)'^^w^+r 


ip'+pr'^' ' 



and 



S = So- /"rf^X J] Qf^,")F«^(") + !!^B{u)-Qa{n)\ ^420) 

r) ^ / 



Hence, m„ are the masses of the d-dimensional vector mesons. To calculate them, one 
needs to solve the eigenvalue problem (|4.16D . 



Solutions of ( [4.16| ) behave at large p like tpip) ^ ci + C2P~^'^~'^^ + ■ ■ ■. Perturbations 
with ci 7^ are non-normalizable; thus we set it to zero and use the shooting method to 
determine the masses of the vector mesons (again, setting d = 3 and n — t- ric). The first 
few masses are 

m^/jL^ ^ 3.08, 15.12, 34.87, 62.32, 97.46, 140.31, 190.86, 249.11, ■ ■ ■ . (4.21) 

They are well described by 

ml/p'^ ^ 3.85^2 + 0.497n - 1.27. (4.22) 

We see that the spectrum of vector mesons is similar to that of the scalar mesons ( [4.11|) , 
except for the absence of an anomalously light meson, as one would have expected. 
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4-3. Pions 

In the D?>/D1 system of section 2, the dynamics breaks a discrete Z2 symmetry 
x^ — 7- —x^. More generaUy, the theory has an 0(6— n) global symmetry, that is dynamically 
broken to 0(5 — n). Unlike the breaking of conformal symmetry, the breaking of 0(6 — n) 
is spontaneous. Hence, one expects to find in the spectrum massles Nambu-Goldstone 
bosons parametrizing the coset 0(6 — n)/0{5 — n); we will refer to them as pions. To 
identify them, recall that in the geometry ( p.l4| ), the Dp-hrane is located at a point on the 
5 — n dimensional sphere S^~^ labeled by fl^-n- The residual 0(5 — n) symmetry is the 
subgroup of the isometry group of the sphere which preserves that point. 

The NG bosons correspond to slow motion of this point on the sphere. Focusing on 
one of them, we can parametrize the AdS^ x S^ geometry as follows: 

ds^ = r'^dx^dx^ + r-^{dp + p^dnl_^ + df + fd9^ + (/ sm9)^dnl_J. (4.23) 

The angle 9 parametrizes one of the directions on the sphere. The NG boson associated 
with it is obtained by promoting ^ to be a function of p, t. Then the DBI action of the 
brane generalizes to 



Sdp = -Id'xjdp-^^^il + p + /20'2)^4 _ ^4J2(1 + J,2)^2. (4.34) 

The quadratic action for 9 is 

^2 = I d'^x I dp[L^E{p)9^ - D{p)9% (4.25) 



where 



Pp-HP + p' 



d— n 
2 



2v/rT7^ (4.26) 

1 



E{p) = lfp--\f + p')'^V^Tr- 

As before, we make a plane wave ansatz 9{p, t) = 9{p)e^'^^^. Then the eom reads 

^^{D{p)9'{p)) + L^mlEip)9ip) = 0. (4.27) 

The general solution of (|4.27|) behaves at large p like a In p + 6. Normalizable modes have 
a = 0. Thus, the mode 9 = constant, which solves ( [4.27|) with mjr = 0, gives a normalizable 
perturbation. 
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Fig. 4: m\l[i^ as a function of fuv^ftvl [x^^'^ in the BKT limit for d = 3. 

As discussed in Section 3, we can separate the threebranes and sevenbranes by chang- 
ing the UV boundary condition to /(A) = fijy, while 6'(A) = 0; this breaks the 0(6 — n) 
rotational symmetry explicitly, and gives the NG boson a mass. The scaling of the pion 
mass with symmetry breaking parameters in QCD is given by the Gell-Mann-Oakes-Renner 
formula 

Flml = a\(^^)\, (4.28) 

where a is the explicit symmetry breaking mass parameter and -F^r the pion decay constant. 
The analog of this in our problem can be read off fig. 4, and takes the form 



m^ ~ 



5.165LVc/y/U - 



d-1 

A\ — 



(4.29) 



As discussed around eq. ( |3.30| ), we can identify a with fuv^ ^ ■ Eq. (|4.29| ) takes then 
the form ( [4.28| ), with ('(/''i/') ~ /x2, Ftt ~ At^~- 



5. Finite temperature and chemical potential 

In this section, we will discuss the effect of turning on finite temperature T and chem- 
ical potential ^ch for the fermion number. We will work at strong coupling and in the 
BKT limit described in section 3. As usual, to construct the phase diagram one needs 
to evaluate the DBI action for various types of D-brane embeddings. We will see that 
at sufficiently high temperature and chemical potential, the symmetric phase is thermo- 
dynamically preferred. Hence, in the Hch ~ T plane there is a line of first order phase 
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transitions which separates the phase with dynamicaUy generated mass and broken sym- 
metry from the symmetric phase. 

The characteristic scale of the temperature of this phase transition is set by the meson 
mass, JI (defined in ( [4.1| )), while the scale of the chemical potential is set by the fermion 
mass //. It will be convenient to rescale the variables 

r^^, /^^, P^j^ (5.1) 

i.e. There = r / L'^ , etc. As a result, the dynamically generated scale is expressed as 

JI=f{0)-Auvexp(-^y (5.2) 



5.1. Finite temperature 

To study finite temperature physics at strong coupling, we replace AdS^ x S^ by the 
AdS-Schwarzschild geometry 

ds"^ = L^ (r'^(F(r)dt^ + dx^dx"^) + !^L, , + dQl] , (5.3) 

Y r^F[r) J 

where a = 1, 2, 3, r^ = vrT, and 

F(r) = l-(^)'. (5.4) 

The induced metric on the probe Dp-brane (|3.1|) is 

dsl, = L^ (r\Fir)dt^ + dxj) + 1 ["^^{p + ff'f + (1 + f")] dp' + ^ciO^i) , 

(5.5) 
where i = 1,2,- ■ ■ ,d — 1. The free energy of the system is given by the DBI action. 



Sc^ Jdp r''--p^-^\JF{p){l + n + ^^(p + ff'Y- (5.6) 

The asymptotic form of the action and equations of motion at energies large compared 
to the temperature {i.e. p ^ rh) are given by eqs. ( p.l6| ), (|3.17| ). The full equations of 



motion are complicated but can be solved numerically. As before, we introduce the UV 
cutoff Kjjv and impose the Dirichlet boundary conditions f{Kuv) = 0. 

In the presence of the black hole, the probe brane can take one of a few qualitatively 
different shapes. One is a "Minkowski embedding" in which the brane does not intersect 
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the horizon. Regularity of the brane at p = leads in this case to the boundary conditions 
( pj.33| ) . The qualitative form of the resulting solution is expected to be the following. At 



low temperatures {vh <^ Jl) , the effect of the black hole on the solution of section 3 is small, 
and in particular the system exhibits dynamical symmetry breaking. As the temperature 
increases, the black hole both grows in size and gravitates more strongly. Eventually, one 
reaches a critical temperature Tc at which the Minkowski solution approaches the horizon. 
As the temperature increases further, the intersection point of the brane and the horizon 
may move. 

In addition to the two types of solution described above, there is always a trivial solu- 
tion, described by /(p) = 0, which corresponds to the phase where mass is not generated 
dynamically and parity is not broken. At low temperatures, we expect the symmetric so- 
lution to have higher free energy than the Minkowski one, and thus be ther mo dynamically 
disfavored. However, above a certain critical temperature, it may have the lowest free 
energy and correspond to the ground state of the system. 

As discussed in section 3, the zero temperature problem is defined by fixing the UV 
cutoff Auv and Kqo (|3.20| ). In the BKT limit one can equivalently fix /I = /(O) ( p^) and 



the ratio C1/C2 in the asymptotic solution 

Turning on finite temperature modifies the full solution, but since the UV boundary 
condition remains the same, the asymptotic form of the solution must still take the formH 

/(p)~p-^(ciln^ + C2y (5.7) 

with the same ratio of C1/C2. 

It is convenient to measure energies in units of the horizon radius rh = tcT, and rewrite 
eq. ( |3.26| ) in the form 



/(p)U=o ~r, ( ^ j ' ^C^ In (^-^ j +C^ln|^^]+ C, 



The constants €> defined in (pT8|) satisfy the constraint 



(5.8) 



— = exp 

rh 



C2 cP' 



(5.9) 



For the rest of this section we specialize to d = 3. 
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To study the finite temperature problem numerically, one can proceed as follows. Integrate 
the full equations of motion, starting with the initial conditions /(O) = JIrp, /'(O) = 0. The 
resulting solution asymptotes to ( pT8| ) for p ^ r^. One can then read off C| ,6*2 and 



use (|5.9|) to calculate p/rh. Finally, one can invert the relation, and determine the jj.rp 
corresponding to the JI of the zero temperature problem. 




Fig. 5: /i/r/i as a function of /i^/rh- 

The resulting graph of Jl/vh as a function of Jlj'/rh is shown in fig. 5. At large 
JLj'/rh, it asymptotes to a straight line with a unit slope, as one would expect, since in this 
regime the effect of temperature is negligible and Jl ~ Jlrp. As the temperature increases, 
one encounters some non-trivial effects. One is that Jl/rh, is bounded from below. This 
implies that the Minkowski embedding does not exists above a certain critical temperature. 
Another is that there is a small region in which for a given Ji/rh there are two distinct 
solutions with different values of Jirp and different free energies. 

We now move on to embeddings of the probe brane which intersect the black hole 
horizon. Such embeddings can be parameterized by the position of the intersection in 



the (p, x^) plane, (po, \/f'h ~ Po)- '^^^ relation between /l/^/i ^ind p^/rh can again be 
determined numerically. The results, which are plotted in fig. 6, exhibit the following 
properties. As po^^O, the black hole embedding asymptotes to the limiting Minkowski 
embedding of figure 5, with Pt = t^h- As po^rh, the solution approaches the symmetric 
one, x^{p) = 0. It is also clear from figure 6 that symmetry breaking black hole embeddings 
only exist in a finite interval of temperatures, rn G (1.51p, 2.43/1). 
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Fig. 6: /i/r/i as a function of po/rh- 

To determine the phase structure of the theory, we need to compare the free energies 
) of the different solutions described above. As mentioned above, these free energies 
are UV divergent, but the differences are finite. It is convenient to subtract from the free 
energy that corresponding to the symmetric solution f{p) = 0. The resulting free energies 
are shown in fig. 7. 




Fig. 7: 5S as a function of th/ p in the BKT limit. The blue curve corresponds to 
the Minkowski embeddings and the pink curve to the black hole embeddings. 



The blue curve corresponds to the Minkowski solution. At low temperatures, this 
(unique) solution is thermodynamically preferred, but for T > Tc ^ O.SG/J/tt, there is a 
first order phase transition to the symmetric solution, /(O) = 0. As the temperature is 
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increased above T = Ti ^ 0.657l/7r, a second Minkowski solution appears, which merges 
with the first one at T = T2 ~ 0.71]I/'k. Only the symmetric solution exists for T > T2. In 
the interval T G {T^,Ti) where T* ^ OA2JI/7:, there exists another type of solutions, which 
originates from the black hole. It is never thermodynamically preferred and represents 
a saddle point of the free energy. It merges with the unstable branch of the Minkowski 
embeddings at T = Ti. 

5.2. Finite chemical potential 

Let us consider the situation where the temperature vanishes, but there is a non- 
vanishing chemical potential conjugate to the fermion number. In the dual description, 
the chemical potential corresponds to the asymptotic value of the temporal component of 
the gauge field on the probe brane. Turning it on gives rise to a modification of the DBI 
action: 

S = jdp r'-y^-'^l + r-^^A',^. (5.10) 

The equation of motion for Aq = ^^^Aq can be integrated to 

"^ ^ ^° -i\ (5.11) 



l + r-A'o 



The conserved quantity d is proportional to the charge density in the boundary theory. In 
order to have a non- vanishing d (and non- vanishing fiux on the worldvolume) , the probe 
brane must intersect the horizon at r = which serves as a source of electric fiux on the 
brane. An important example of a solution with fiux is given by the straight embedding 
/(p) = 0. Of course, it is still possible to have a Minkowski embedding where the D-hrane 
does not intersect the horizon. In this case the electric field on the worldvolume vanishes 
and Aq = Hch everywhere. 

One can write Hch = Ao{p^oo) as an integral 

7-2 T-2 i-oo 

2ti 2n Jq 

where the boundary condition Aq^t = 0) = is ensured by requiring the temporal Wilson 
loop to be well defined at the point where the Euclidean time circle shrinks to zero size. 
This also ensures that at vanishing Hch the charge density vanishes. One can now relate 
d to the chemical potential fich- For the straight D-brane embedding, f{p) = 0, one can 
solve for Aq using ( |5.11| ) and then integrate ( |5.12| ) to obtain 'JlJ^~ = {4r {5 / 4)'^ /^/n)d. 
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Fig. 8: 5S as a function of /I^^Z/ig: straight pink line - Minkowski embedding with 
^0 = IJ-ch] blue curve - f{p) = embedding; yellow curve - solution with flux which 
interpolates between the two. 

To obtain the phase diagram, let us first consider the Minkowski embedding and the 
straight /(p) = embedding with flux. We will fix /i = /(O) for the Minskowski embedding 
and will measure energies in units oijl (note that according to (|5.12|) , the natural scale for 
lich is /i = vA/i, the scale of the fermion mass). The action of the Minkowski embedding 
is constant, Sm = — 0.0478/i'^. For the straight D-brane embedding, one can use ( p.lO|) to 
show that S'/=o = —0.09771^/^. Clearly, there is a first order phase transition between the 
Minkowski phase with vanishing density and the phase where the D-brane is straight and 
has a finite density. The complete phase diagram, shown in fig. 8, is a little more involved, 
due to the existence of a third type of solutions which have non-vanishing worldvolume 
flux and intersect the horizon, but have a nontrivial profile /(p); they are discussed in 
Appendix B. 

We have not analyzed the phase diagram with both temperature and chemical poten- 
tial non-vanishing. However, as in other D-brane systems, we expect a smooth line of first 
order phase transitions in the Hch — T plane. This line presumably intersects the p = 
line at Tc ~ /I and the T = line at fj.ch,c ~ P, as discussed above. 



6. QCD 

Four dimensional QCD is expected to undergo a conformal phase transition as one 
varies the number of colors and fiavors. The dynamics in the vicinity of this transition may 
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lead to a phenomenologically viable model of dynamical electroweak symmetry breaking, 
known as "walking technicolor" (see e.g. Jl^ for a review). Since the transition occurs in 
a regime where the theory is strongly coupled, it is difficult to analyze its dynamics there 
using conventional field theoretic methods. One may hope that holographic techniques 
may shed light on it. In this section, we briefiy review the expected structure near the 
phase transition from the point of view of our analysis. 

We start with a gauge theory with gauge group G = SU{N) and F fiavors of Dirac 
fermions in the fundamental representation of the gauge group (see e.g. [^ for a review). 
This theory is infrared free for F > UN/ 2; for smaller F it is asymptotically free in the 
UV. As we review next, its IR dynamics depends on the number of fiavors and colors. 

If F is only slightly smaller than llA^/2, the gauge coupling runs from zero in the UV 
to a small non-zero value in the IR |^1[]. The theory dynamically generates a scale, Aqcd, 
which can be thought of as the crossover scale between the UV and IR. For energies well 
above this scale, the correlation functions are dominated by their (asymptotically free) UV 
forms, while for E ^ Aqcd they are dominated by the interacting IR fixed point. 

It is important to stress that the IR fixed point is an isolated CFT, which does 
not depend on any continuous parameters, unlike the defect theories studied in previous 
sections, which depended on the continuous parameter A, the 't Hooft coupling of A/" = 4 
SYM. Instead, the QCD fixed point depends on the discrete parameter x = N/F. We will 
treat this parameter as continuous by taking A^ and F to be large. 

As F decreases (or x increases) , the IR theory becomes more strongly interacting. In 
particular, the dimensions of operators deviate more and more from their free values. One 
can think of x as a measure of the strength of the interactions in the IR CFT. Thus, it 
plays a role similar to that of A in the defect CFT discussed in previous sections. 

When the coupling exceeds a critical value x = Xc, the model undergoes a phase 
transition. For x > Xc, some of the symmetry of the IR CFT is dynamically broken. In 
particular, the SU{F)l x SU{F)fi global symmetry, which is a symmetry of the Lagrangian 
preserved by the quantum dynamics for x < Xc, is spontaneously broken to the diagonal 
SU{F). The conformal symmetry of the IR CFT is broken as well, at the scale /x ~ 



AQCDexp(— a/i/x — Xc) mentioned in the introduction (see fig. 9). 
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Fig. 9: Conformal symmetry breaking scale /i as a function of the number of 
flavors F. 

This transition is expected to be of BKT type. One can probe it by studying the vac- 



uum expectation value of the operator OW ( p.4| ). For x < Xc, the IR theory is conformal, 
and this vev takes the form ( p.5|) . For x > Xc, the IR conformal symmetry is violated by 
a mechanism similar to that of section 2. As x — ?■ Xc, the anomalous dimension 7 of i/ji/j 
approaches —{d — 2)/2 = —1 [jl^, and beyond that point it formally becomes complex. The 
interpretation is the same as in the above discussion of the defect theory. At x = Xc, two 
fixed points of the RG which differ by the coefficient of the quartic double trace operator 
[t/jt/j)'^ merge and annihilate (or move off to the complex plane). This generally leads to 
Miransky scaling of the dynamically generated scale [|12| . 

In the analysis of the earlier sections, we had a conformal field theory for A < Ac. The 
analog of that here is the infrared fixed point of QCD for x < Xc- The RG fiow of the gauge 
theory from its free UV fixed point to that IR theory corresponds from this point of view 
to a particular choice of UV cutoff. Hence Aqcd corresponds in our previous discussion 
to the UV cutoff Ajjv (|2.31|) . The strength of the interactions of the fermions at the fixed 
point, which in our previous discussion was controlled by the 't Hooft coupling A, is in 
QCD controlled by the parameter x [i.e. the number of colors and fiavors) rather than the 
(running) coupling. To focus on the physics of the phase transition one can consider the 
double scaling limit x — )■ Xc, Aqcd -^ 00, with the dynamically generated scale fi held 
fixed. This is the analog of the BKT limit discussed in previous sections. 

It is interesting to ask what our analysis of the CPT in defect theories teaches us 
about the corresponding phase transition in QCD. One question concerns the nature of 
the second fixed point that merges with the IR fixed point of QCD at x = Xc. For x slightly 
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below Xc, this fixed point can in principle be obtained by adding the (slightly irrelevant) 
double trace operator {ifjil^)'^ to the Lagrangian, and following the resulting RG trajectory 
to its UV fixed point. Of course, in practice this is hard since the theory one perturbs is 
strongly coupled. 

As X decreases, the distance in coupling space between the two fixed points increases. 
An interesting question is whether the second fixed point exists all the way to weak cou- 
pling. In our analysis of defect theories, this was the case: setting A = and adding 
the irrelevant four Fermi perturbation led to a non-trivial UV fixed point at which the 
dimension of the fermion bilinear was equal to one (|3.9|) , in agreement with |3^ . 

Thus, it is natural to conjecture that in four dimensional QCD the second fixed point 
also persists to arbitrarily weak coupling. There is a small subtlety here: the dimension of 
the fermion bilinear at the second fixed point is given hj d — 1 + 'j- ( |3.9| ) . As the gauge 
theory coupling approaches zero, which in QCD happens when x — ?■ 2/11, this dimension 
approaches one. So at x = 2/11, the operator i/ji/j becomes a decoupled free field, and 
adding (i/ji/;)'^ to the Lagrangian simply gives mass to this field. In particular, it does not 
take one back to the free field theory one started with. This problem is related to the fact 
that the NJL model is not renormalizabe in (i = 4, but rather has a logarithmic residual 
dependence on the UV cutoff. 

The above difficulty is avoided for any non-zero coupling, since then the analog of 
C(A) in (|3.9| ) is strictly positive and the dimension of the fermion bilinear is larger than 
one. Thus, in QCD it is natural to conjecture that there is a new fixed point obtained by 
adding the four Fermi perturbation to the Lagrangian for any non-zero value of the gauge 
coupling. 

In the defect problem, that fixed point can be analyzed via standard field theoretic 
techniques, described e.g. in |^^. It is interesting to ask whether one can establish its 
existence in four dimensional QCD as well. There are two complication in doing this. One 
is that as we just mentioned, one needs to turn on a non-zero gauge coupling. This is not 
a major problem since the fixed point should exist for arbitrarily small gauge coupling, 
and therefore should be accessible to standard weakly coupled field theory techniques. 

A more serious issue is that the fermions "0 in QCD are large A^ matrices, while in 
the defect problem they are vectors. More precisely, in both cases they are matrices of size 
N X F, but while in the defect problem the transition takes place when A^ is large while F 
is of order one, in QCD it happens at a finite value of x = N/F. As is well known, while 
vector models are typically solvable at large N using the methods reviewed in [^2|], matrix 

39 



models are in general much more difficult to solve. Thus, the study of weakly coupled 
gauge theory with a four Fermi interaction is an open problem. It is natural to expect 
that this model does exhibit a non-trivial UV fixed point in which the dimension of the 
fermion bilinear is given by (|3.9|) . 

A widely discussed (and controversial) issue in the context of the CPT in QCD con- 
cerns the mass of the dilaton. The question is whether in the BKT limit, the mass of 
the lightest massive scalar goes to zero or remains comparable to the masses of the other 



mesons. This is an old debate; some recent contributions appeared in [^ , |36[] . 



It turns out that this issue is directly related to the vev of the trace of the stress 



tensor T^. The authors of |^ used the fact that (T^) ~ /3(/i)((9), where /3(/i) is the beta 



A* ' 



function for the coupling h of an operator (9, which in this case is a four fermi operator. 
They argued that the natural size of the vev of O is /U^, which leads to (T^) ~ /3(/i)/U^. 
Since the beta function goes to zero in the BKT limit x — )■ Xc they concluded that the 
expectation vaue of T is parametrically suppressed relative to the natural mass scale of 
the theory, fi. One can then show that the mass of the dilaton is suppressed as well. 
The authors of [^ argued that in fact the expectation value of the stress tensor takes its 
natural size, (T^) ~ //^, and the mass of the dilaton is correspondingly not suppressed. 

We calculated both the expectation value of the stress tensor and mass of the low 
lying a mesons for the defect theories. We found that the stress tensor (|3.31| ) and mass 



of the lightest scalar meson ( [4.1C1| ) are not parametrically suppressed, in agreement with 
3^ . As we mentioned in section 4, the latter does exhibit a mild finite suppression, which 



may also occur in the BKT regime of QCD and play a role in walking technicolor. In 



particular, the authors of [^ found that the mass of the lightest scalar is smaller than 
that of the lightest vector meson by a factor of about 2.8. In our system this ratio is (see 
( [4.10|) , (|4.21|) ) about 2.6. Of course, these numbers cannot be compared directly, since the 
analysis of [^ is not done in a controlled approximation, and ours is done in a different 
system. However, it is intriguing that the two calculations yield similar suppressions. 

Finally, we would like to emphasize that the conformal nature of the phase transition 
at a; = Xc in QCD is only a scenario, which has not been proven to date. It assumes 
that the transition has a sensible limit as Aqcd -^ oo? ^-c it is continuous. An alternative 
scenario is that the system undergoes a first order phase transition in which the dynamically 
generated scale jumps from zero to a scale of order Aqcd at a value of x for which the 
dimension of '^'0 is larger than two. Of course, in that case the discussion of this section 
is inapplicable. 
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7. Discussion 

In this paper, we showed that the system of defect fermions coupled to A/" = 4 SYM 
( p.2|) undergoes a phase transition between a weak coupling phase, in which the theory is 
scale invariant, and a phase in which the fermions acquire a dynamical mass. The order 
parameter exhibits Miransky (or BKT) scaling (|2.25|) , a behavior also expected to hold in 
QCD near a critical number of flavors. 

At flrst sight, it is surprising that this type of continuous, inflnite order transition 
is stable to large changes in the parameters of the theory. In other examples of systems 
with holographic duals, even when a transition is continuous at weak coupling, it typically 
becomes strongly flrst order when one goes to strong coupling. It seems that what pro- 
tects the conformal phase transition in our system is its topological nature. It provides 



a realization of the picture proposed in |12], according to which CPT's occur when two 
flxed points of the RG approach each other and annihilate as one changes the parameters 
of the model. This mechanism is robust - as long as the physics varies continuously with 
the parameters, a transition that is conformal in one regime will remain so. 

In the class of systems ( |3.3| ), we saw this mechanism in action. By changing the 
parameters d and n, we were able to push the transition to weak or strong coupling, but 
since the physics is expected to vary smoothly with these parameters, it is natural that it 
was found to be conformal in both regimes. 

Since the theory develops an arbitrarily small mass scale near the transition, it is 
interesting to ask whether the physics in this regime is universal. We provided descriptions 
of that physics at weak and strong coupling. In fleld theory, we discussed the order 
parameter M{p) (|2T0| ), ( pT5| ); its analog in gravity was f{p) (g]29|), QOOD . The two 
order parameters satisfy very similar second order linear differential equations ( p.4|) , ( |3.18|) , 
and encode the properties of the fermion bilinear near the transition at weak and strong 
coupling. An important difference between the two is that while in the strong coupling 
(gravity) regime we provided a precise deflnition of /(p), as a scalar fleld describing the 
position of a probe brane in the bulk, in the fleld theory analysis we only deflned M{p) 
as the fleld measuring the deviation from conformality to leading order, ( |2.1C1| ). A natural 
question is whether one can deflne M{p) beyond the leading order and for all values of the 
coupling, such that for strong coupling it approaches /(p). 

As a step in this direction, we would like to offer the following observation. Consider 
the action 

7d_ j„ P \ „d-l 



Snp = I d'x I dp( -^^ ] /- Vl + f"- (7.1) 

41 



This action was obtained by replacing n—d in ( |3.16| ) by C(A) (and omitting an unimportant 
overall constant). The relation (|3.19|) guarantees that (|7.1j ) gives the correct action for all 
d and n in the strong coupling limit A — ?■ oo. However, there is no a priori reason for it to 
be reliable away from strong coupling. 

To see how well it does, it is interesting to take the opposite limit, A — > 0, where 
C(A) is small (|3.5|) and one can expand in it. The equation of motion of ( |7.1| ) in that limit 
can be read off ( |3.17| ) with the replacement n — d ^ C. To compare to the perturbative 
analysis of appendix A, we take /(p) = m + 0(A), and omit terms of order A^ or higher. 
This leads to: 

|(/-/VC(A)A'^^;^^a (7.2, 

Remarkably, this agrees precisely with the weak coupling result (A. 15), which was obtained 
by evaluating one loop Feynman diagrams! If this agreement is not accidental, it might 
indicate that the action ( [7.1D provides a description of the physics near the phase transition 
for all values of d and n. In particular, it may be valid for the original D3 — D7 system 
discussed in section 2, when C(A) is close to 1/4. 

Of course, the action (|7.2|) depends on C(A), which is related to the anomalous di- 
mension of t/jt/j, and is only known in the two limits A — ?■ 0, oo. However, its precise form 
is less important than the fact that it passes through the critical value [d — 2)^/4. The 
physics near the transition depends on k (|3.6|) , and ultimately on fx ( |3.25|) . 

Although the discussion above is suggestive, there is clearly much to understand re- 
garding universality in this class of systems. In particular, it would be interesting to show 
that when the dynamically generated mass is small relative to the UV cutoff, the physics 
of the order parameter is given by the action ( [7.1| ) without a' corrections. It would also be 
interesting to include in the discussion the gauge fields dual to flavor currents, and extend 
the definition of the order parameter M(p) in field theory beyond leading order. 

Other issues worth exploring include: 
(1) We focused on the case of one flavor brane (F = 1). One qualitatively new effect 
that appears in the D3 — D7 system for F > 1 is the breaking of the U{F) flavor 
symmetry to U{Fi) x U{F2) with Fi +F2 = F. In the brane picture this corresponds 
to a state in which Fi branes are deformed to x^ > and F2 to x^ < 0. The 
spontaneous breaking of the symmetry leads to the appearance of Nambu-Goldstone 
bosons which parametrize the coset U{F)/U{Fi) x U{F2). They can be studied as in 
|T8[] , by analyzing the relevant components of the gauge field Ap. 
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(2) The model we studied is closely related to A/" = 4 SYM, about which much is known 
due to its integrability at large A^ (see e.g. [^2| for a review). The defect fermions 
can be viewed as replacing closed spin chains by open ones, with particular non- 
supersymmetric boundary conditions. It might be possible to use the description in 
terms of open spin chains to say more about the model, and in particular determine 



the function C(A), and perhaps understand the role of the action ( [7.1|) for general A. 

(3) Our discussion took place at large A^, and it might be interesting to consider 1/A^ 
corrections. Beyond the leading order many aspects of our discussion need to be 
reconsidered. In particular, it is not obvious that the model is still conformal even at 
arbitrarily small A. The A/" = 4 SYM coupling in the bulk of course does not run, but 
one can in principle have a non-zero /3-function for the boundary coupling in ( |3.3| ) due 
to loops of fermions. 

(4) It might be interesting to generalize the discussion to other backgrounds, such as linear 
dilaton spacetimes. In these spacetimes there is an analog of the BF bound, and they 



exhibit holography |^,Q. In this context it is useful to mention two dimensional 



string theory (or c = 1 matrix model) [^ and the 3 + 1 dimensional LST of the 
conifold (see e.g. |^3) which involve condensates of closed string fields living at the 
boundary between the stable and unstable regions. 

(5) In |I[] it was suggested that the D?> — D7 system ( ^.1| ) might provide a model for the 
interactions of electrons in graphene. It would be interesting to realize the conformal 
phase transition of this system in this or other condensed matter system. 

(6) It would also be interesting to extend the picture presented in this paper for interacting 
defect fermions to QCD near the critical number of flavors Fc. 
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Appendix A. Weak Coupling Analysis 

In this appendix we calculate the one loop anomalous dimension 7(A) of t/jt/j in the 
gauge theory (^j^) , and derive a nonlinear generalization of the differential equation 
for M{p). 

Adding a fermion mass term to (|3.3|) leads to the action 



S = Sn=4 + / d'^x{iip^i; + gi;^^; + gtl)TM(t)^^ - mi^i^). (A.l) 

Here and below, the indices (a, b) and (M, A^) run over the d directions along the defect 
and (6 — n) transverse dimensions respectively. The Feynman rules on the defect are: 

l^ermion propagator : 



q2 _ j^2 



Q , ^7-10 I d q± iSmn ■ Smn -'^("2" 

Scalar propagator : D^ — ' — . 






Gauge field propagator : Dl^^ = - 2^/2 ^^^ 9ab - {d - 2)g„gb(l - 0) ^4_d^^^_d)/2 ^ 

vertex : igjat\ igTMt\ 

(A.2) 

where the scalar and gauge propagators are obtained by integrating the four dimensional 

propagators over the (4 — d) dimensions transverse to the intersection, q±; t* are SU{N) 

matrices. As discussed in Section 3, Tm and 7^ can be thought of as SO{d + 5 — n, 1) Dirac 

matrices. In particular, {r,7} = 0, and in mostly minus signature {r^''^,r^} = 2r]^'^^ = 

The one-loop fermion self energy 



^S(p) 



r 


■^ 


1 


1 




(A.3) 



is given by 

-zE(p) = 



A /• d'^k 



2 } (27r) 



'^^^(^-^)^^ra^^+^"'^(^-^)^"ra^^' 



(A.4) 
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Using the parametrization of E in (|2.13| ), one has: 



„, , iX f d'^k 



DlAp-k)l' 



m 



n 



m 



■*w - ^ + 1''/ (S-? K'^ - '^■''w^-'' ^ ^-(^ - *'r°ra^' 






After Wick rotating and inserting the propagator, one finds 



(A.5) 



A r(^) /• „, 1 1 



1 



^^^ 2 167r(4+rf)/2 J " "^2 |p _ ^|d-2 )^2 ^ ^2 



(A.6) 



1 

X tr- 
d 



((^ - 2) + (6 - n))^^ + 2{p-kY 



the second of which can be rewritten as 

Aiv) - 1 + ^^^S^ /rf^fcM^^l±i^^ 
^^> " ^ 2 167r(4+'^)/2 J ""^ 2(A;2 + m^) 



k' 



+ 



{d-2){l-0 
4(A;2 + m2) 



|p_/j.|d 2 p2|p_/j.|d 2 p2|p_/j.| 

'.4 .2 2A;2 ^ 



d-4 



|p — /c 



|d-2 



p/|p 



k^ k^ p2 2A;2 

1 \ 

_ ^|d-2 p2|p_^|d |p _ /j|d |p_/j|d 



Now. 



recaU the identity for the Green's function of the Laplacian: 



J___(rf-2)27r^ (,) 

d-9 -n/r/N " V 



V" - 

j,d—2 Pf' — ^ 



r). 



(A.7) 



(Ai 



Applying it to (p^), ( [A77| ) yields 



9p 



|,(/-'B'w) + li^ O' 



?7l 



— 1-0 + 6-n U"-^ ^ =0, 

2 / p + m^ 



„d-l 



|(/-'^'(rt) + (, - 2) ((6 - n -Ki - 2)gA_ _ (l_^,l _ Og;j^A_J p»-.__ 

(A.9) 
To leading order in m, eqs. (|A.9| ) are solved by 



= 0. 



B{j)) =m(l + /31n| — j j ~m| — 






p 



(A.IO) 
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where the coefficients are 



„ = _Ar(a-n + .-2)i-il^(i-0. ^ 



87r2 V 'd 2 ^ ^'d-2^ 

and /i is the renormahzation scale. Note that in the absence of the scalar fields [i.e. setting 
71 = 6) and for (i = 4, ct is proportional to the gauge parameter ^, in agreement with |J2^. 
As another check, the one-loop values of a, (3 in four-dimensional Yukawa theory and QED 
are well-known results, derived for example in [i^]. Our results agree with theirs under 
the map g'^ o A/ 2. 

The quantities derived in ([A.1C1| ) can be combined to form the gauge-invariant part of 
the fermion self-energy ( |2.15| ): 

where 

7(A) = /3 - a = -^^^^4^ (2(d - 1) + (6 - n)) . (A.13) 

7(A) is the anomalous dimension oiifji/j, (|2.7|) . It is gauge invariant, as expected. 

In section 3, we saw that M{p) satisfies eq. ( p.4|) at large momentum. That equation, 
and its nonlinear generalization which extends to all p, can be derived to order A from ( [A.9| ) 
in the following way. Expanding A{p) ~ 1 + Ai{p) + O(A^), B{p) ~ m(l + Bi{p) -\- 0(A^)) 
and using ([A.12|) , one finds 



M{p) ~ m(l + Bi{p) - Ai{p)) + O(A^). (A.14) 



If we now take the difference of the equations in (|A.9|) , we arrive at 



whereM 



|-(/-iM'(p)) + C(A)/-i^^ + 0(A2) = 0, (A.15) 



C{\) = -(d - 2)7(A) = 4^(2(ci - 1) + (6 - n)) + 0{X'). (A.16) 
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Note that this is consistent with eq. ( |3.8| ) when C(A) is small: there we can expand 



2 \d-2j ' W d-2 
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In particular, the {d,n) = (2,6) result for C(A) agrees with [jT2|, and the {d,n) = (3,5) 
result for the D3/D7 system of section 2 is 



C{X) 



5A 

12^ 



+ 0{X' 



Finally, we replace m with M{p), which is consistent at leading order in A, and find: 

-1 Mip) 



d_ 
dp 



{p'-'M'{p)) + C{X)p'^- \,^'M[^y + 0{X') = 0. 



(A.17) 



For large p (or small M) this reduces to (|3.4| ) . 

Eq. ( |A.17| ) describes the vicinity of the conformal phase transition at weak coupling. 
It is an improvement over ( |3.4| ) which breaks down at momenta on the order of the dy- 
namically generated scale, and can be thought of as the analogue of the nonlinear strong 
coupling equation (|3.17| ) following from the DBI action. 

Near d = 2, where weak coupling techniques are applicable, ( |A.17|) can be alternatively 
obtained by using the Schwinger-Dyson equations in the rainbow approximation. These 
are given diagrammatically by 



iE(p) = 



O^ 



+ 




(A.18) 



where the lines with grey dots are full (as opposed to bare) propagators, and the black 
dots denote the fuU vertices A^ = igV" + 0{g^), A^ = ig'yf' + 0{g^). Then Eq. (^) is 
promoted to 



-^np) = -\j 



X r d'^k 



{2ny 



DUp-k)r 



A{kW + B{k) ^, , .^^ ,^,^„ A{k)l^ + B{k) ^, 



A{kYk'^-B{kY 



K + Dl,{p-k)V^ 



A{kYk^-B{kY ^ 
(A.19) 



which, again, can be separated into different spinor structures: 



, , iX f 
B[pj = m tr / 



d'^k 
(2^ 
^. ^ , ^A [ d'^k ^ 



DZAp-k)Y 



B{k) 



DUp-kh 



A{kYk^-B{kY ^ 



K+Dt,{p-k)T'^ 



B{k) 



A{kYk^-B{kY 



A'; + d:,{p - k)^^ 



A{kyk^-B{kY 

A^ 



A{kYk'^-B{kY 
(A.20) 



At 
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Now by eq. ( |A.11| ), the deviation of A{p) from unity is suppressed by {d — 2) relative to 



B{p). It is therefore consistent near d = 2 to set A{p) = 1, and replace B{p) with M(p). 
Likewise, the gauge-dependent part of the gauge field propagator ( [A .21 ) is suppressed by 
(ci — 2), so we are free to fix ^ = 1 in the following. Finally, to order A, it is a consistent 
truncation to work in the quenched approximation and replace the full vertices with bare 
ones. 

With these approximations and after Wick rotating, the first line of ( |A.20 ) becomes 



Upon applying the Laplacian, we arrive at the differential equation 

The deviation of ( |A.22| ) from ( [A.17| ) appears at order {d — 2), and presumably is reconciled 
when one takes higher order corrections into account. 



Appendix B. Some holographic results in the BKT limit 

B.l. Action 



Let us investigate the behavior of the action ( |3.16D in the BKT regime. It is intuitively 



clear that in the limit k— t-O the action approaches a constant whose scale is set by /i. We 
will see momentarily that this is indeed the case. It is technically convenient to divide the 
region of integration in ( |3.16| ) into the small p part and the large p part. Let them be 
separated by p*, such that p <^ p^ <^ ^uv, so that the total action consists of two terms, 

S[f] = S[f]o<p<p. + ^[/]p.<p<A,.. (B.l) 

In the large p region the action is approximately quadratic in /(p), f'{p) and therefore one 
can use equations of motion to show that the non- vanishing contribution to the integral 
only comes from the boundary term at p = p*: 

S[fU<,<A,, = -\ [pTip]]' f{p)\p=p,. (B.2) 

The boundary term coming from the UV cutoff vanishes precisely because /(p = h.uv) = 
0. The integral from p = to p = p* can be evaluated numerically. The result of 
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this numerical integration, 5'[/]o<p<p*7 with the boundary term (|B.2|) added, converges 
to — C/U"^V3, where C is a positive numerical factor and V3 is the volume of the three- 
dimensional spacetime. This implies, according to ( |3.25| ), that 

5o~-A'exp(-^) (B.3) 



consistent with the BKT-type phase transition. Here the zero in the subscript indicates 
the vacuum nature of the solution. 

B.2. Chemical potential 

There are three types of embeddings: Minkowski embeddings with vanishing gauge 
flux on the woldvolume of the D-brane (so that Aq = fXch), straight embeddings f{p) = 
with flux, and curved embeddings with flux which start at p = and can be parameterized 
by the initial condition /'(O). The flrst two were discussed in Section 5.2. 

To analyze embeddings with flux and nontrivial proflle, it is technically convenient to 
measure all energies in the problem in units of d. For a flxed value of d, one can vary /'(O) 
to obtain the solutions with various asymptotics (in practice, one integrates the equation of 
motion for /(p) which follows from (|5.10|) while eq. ( |5.11|) is used to eliminate Aq in favor 



of d). The value of Jlch/d can then be obtained by integrating ( |5.11[ ) numerically. Finally, 
the value oiji/d for a given solution can also be obtained from the asymptotic behavior, in 
a similar way as the values oiji/rh and JL^/rh were obtained in the flnite temperature case, 
where eq. (|5.9|) was used. One can therefore compute the ratio Jich/J^ ^s a function of /'(O). 
This turns out to be a monotonically increasing function, which starts from 'pch/'P ~ 0-6 
(which corresponds to the /'(O)— t-O limit, where the solution joins the straight / = 
solution with flux) and approaches unity from below (which corresponds to /'(O)— t-oo, the 
regime where the curved solution with flux joins the corresponding Minskowski solution). 
The action can be computed numerically; the results are shown in Fig. 8. These curved 
solutions with flux are never thermodynamically preferred. 
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